Soliton-dynamical approach to a noisy Ginzburg-Landau model 
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We present a dynamical description and analysis of non-equilibrium transitions in the noisy 
Ginzburg-Landau equation based on a canonical phase space formulation. The transition path- 
ways are characterized by nucleation and subsequent propagation of domain walls or solitons. We 
also evaluate the Arrhenius factor in terms of an associated action and find good agreement with 
recent numerical optimization studies. 
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Phenomena far from equilibrium are widespread in- 
cluding turbulence in fluids, interface and growth prob- 
lems, chemical and biological systems, and problems in 
material science and nanophysics. Here the dynamics of 
complex systems driven by weak noise, corresponding to 
rare events, is of particular interest in the context of e.g., 
nucleation during phase transitions, chemical reactions, 
and conformational changes in macromolecules. The 
weak noise limit is associated with a long time scale cor- 
responding to the separation in energy scales of the ther- 
mal energy and the energy barriers between metastable 
states; the transition takes place by sudden jumps be- 
tween metastable states followed by long waiting times 
in the vicinity of the states. The fundamental issue is 
thus the determination of the transition pathways and 
the associated transition rates. 

A particularly interesting non equilibrium problem of 
relevance in the nanophysics of switches is the influence 
of thermal noise on two-level systems with spatial de- 
grees of freedom, see |], ||, 0. In a recent paper by 
E, Ren, and Vanden-Eijden this problem has been 
addressed using the Ginzburg-Landau equation driven 
by thermal noise. These authors implement a powerful 
numerical optimization technique for the determination 
of the space time configuration minimizing the Freidlin- 
Wentzell || action and in this way determine the orbits 
and their associated action yielding the switching proba- 
bilities in the long time-low temperature limit. The pic- 
ture that emerges from this numerical study is that of 
noise-induced nucleation and subsequent propagation of 
domain walls across the sample yielding the switch be- 
tween the two metastable states. 



In recent work we have addressed a related problem in 
nonequilibrium physics, namely the Kardar-Parisi-Zhang 
equation or the equivalent noisy Burgers equation de- 
scribing for example a growing interface in a random 
environment. Using a canonical phase space method de- 
rived from the weak noise limit of the Martin-Siggia-Rosc 
functional Q or directly from the Fokker-Planck equa- 
tion J^, [| , we have in the one dimensional case analyzed 
the resulting coupled field equations minimizing the ac- 
tion. The picture that emerges is that the transition 
probabilities in the weak noise limit are associated with 
soliton propagation and nucleation resulting from soliton 
collisions. 

In this letter we apply a soliton approach in the canon- 
ical phase space formulation to the noisy Ginzburg- 
Landau equation and attempt to account for some of the 
numerical results of E et al. We thus give analytical argu- 
ments for the propagation of noise induced domain walls 
or solitons, the nucleation events associated with domain 
wall creation and annihilation, and the associated time 
dependent action. Details of our analysis will be given 
elsewhere. 

The noisy Ginzburg-Landau equation for a field u(x, i) 
driven by white noise has the form 

S = " r £ + ?? ' fafo*M '°)> = &S(x)6(t), (1) 



with free energy 
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In the switching problem considered by E et al. J| V(u) 
is given by the "Mexican Hat" double well potential 



V(u) 
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with strength parameter kg. T is a kinetic transport co- 
efficient setting the time scale. 

The Ginzburg-Landau equation in its deterministic 
form has been used both in the context of phase order- 
ing kinetics fic[| and in its complex form in the study 
of pattern formation fll|| . In the noisy case for a finite 
system the equation has been studied in fl2|| ; see also an 
analysis of the related </> 4 theory in jl3| . In the present 
problem the noisy equation provides a generalization of 
the classical Kramers problem [O to spatially extended 
systems. 

The equation admits a fluctuation-dissipation the- 
orem yielding the stationary distribution P s t a t oc 
exp[— 2TF/A] . The equilibrium states follow from 
8F/Su = —d 2 u/dx 2 — 2k 2 u(l — u 2 ) — 0, giving the two 
degenerate uniform ground states u = ±1 with F = 0, as 
well as nonuniform domain wall solutions 



Ud-w{x) — ±tanhfco(£ — #o)i 



(4) 



centered at xq, connecting the two ground states. The 
associated free energy is 

F dw = 4fc /3. (5) 

Since the spectrum of 5 2 F/Su 2 = -d 2 /dx 2 - 2fcg(l -3u 2 ) 
is positive for the ground states and has zero-eigenvalue 
Goldstone modes (translation modes) for the domain wall 
solutions, see e.g. we infer that the free energy land- 
scape possesses two global minima at u = ±1 and a series 
of local metastable saddle points of free energy 4nfco/3, 
corresponding to n connected domain walls at positions 

X i ■ % — 1 • * " " T\i • 

In order to address the issue of noise-driven transi- 
tions in the Ginzburg-Landau equation (^) we apply 
the phase space approach developed for the Kardar- 
Parisi-Zhang equation || This method is based on 
a weak noise WKB-like approximation, P({u(x)},i) oc 
exp[— S({u(x)}, t)/A], applied to the Fokker-Planck 
equation AdP/dt = HP for the transition probability 
P({u(x)}, t), driven by the Hamiltonian or Liouvillian 
H({u(x)}, {d/du(x)}). To leading order in A the ac- 
tion or weight function S then satifics the Hamilton- 
Jacobi equation dS/dt + H = with conjugate mo- 
mentum or noise field p(x) = SS/Su(x). Moreover, 
the underlying principle of least action SS = yields 
Hamiltonian equations of motion du/dt = SH/Sp(x) and 
dp/dt = —5H/Su(x). We note that the equations of mo- 
tion are identical to the saddle point equations in the 
Martin-Siggia-Rose functional formulation J(| [5J . 

For the Ginzburg-Landau equation we then obtain ex- 
plicitly the coupled deterministic field equations 



du 
~dt 
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T^ + 2Tk 2 u(l-u 2 )+p 7 



(6) 



f = -rg-2L, o yi-3A 



derived from the Hamiltonian 



H = 



dxp[p + 2L— + ATk 2 u(l - u 2 ) 



(8) 



In the weak noise limit the transition pathways are de- 
termined by the orbits in the (u,p) phase space from an 
initial configuration u\(x) at time t — to a final config- 
uration u^ix) at time t = T. The conjugate field p is a 
slaved variable representing the noise driving the system. 
Finally, the transition rate P(u\ — ► U2,T) oc exp[— S/A] 
is determined by the action, H = j dxTL, 



S(ui 



pu 2 ,T 

U2, T) = / dxdt 
J «i,0 



du 



(9) 



associated with the orbit. 

To linear order, the effect of the noise field p is to im- 
part a velocity to the static domain wall (||) . This is seen 
by expanding u and p on the translation mode associated 
with the static domain wall, u tm = —(l/m)dud v/ /dx, 
where the mass m is to be determined. Setting u = Udw + 
UoUtm and p = p Q u tn i and using (S 2 F/Su 2 ) Udv u tm = we 
obtain from (|^) and (0) duo j dt = po and dpo /dt = with 
solutions uq — pot and po = const. Consequently, 



u(x, t) ~ u dv 



Pot u dw 
m dx 



Udw (x - —t) , (10) 
V ml 



describing a domain wall of mass m propagating with 
velocity v — po/m. With the above normalization of 
the translation mode the noise field po is a momentum 
contributing to the total momentum n = J dxudp/dx 
(generator of translation)and thus canonically conju- 
gate to the position x of the domain wall. From 
J da; cosh - k$x — 4/3fco we infer the mass m. For the 
energy and action associated with the propagation of a 
single domain wall in time T we then have 



4fci 



2m 2m 6 



(11) 



This analysis generalizes directly to a dilute gas of con- 
nected non-overlapping domain walls. The noise p then 
gives rise to individual velocities imparted to each do- 
main wall and the energies, momenta, and actions simply 
add up. 

In addition to the translation modes, time-dependent 
diffusive modes are also excited, corresponding to small 
Gaussian fluctuations about the local minima. In the 
neighborhood of u = ±1 it follows from (^) and (^) that 
u and p in a plane wave decomposition, exp(ikx) evolve 
with time-dependence exp(±r(fc 2 + 4k 2 )t). About the 
local saddle points corresponding to domain wall propa- 
gation, the extended modes are phase shifted correspond- 
ing to the trapping of a localized deformation mode with 
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time dependence exp(±r3fc^t) and the translation mode 
discussed above, see e.g. [[15). 

The dynamical interpretation of the noise-induced 
switching in the Ginzburg-Landau equation is now clear. 
The transition from e.g. u = +1 to it = —1 starts by 
nucleating a small region of size ~ fc^ 1 and then prop- 
agating a domain wall or walls, with superposed linear 
modes subject to energy and momentum conservation 
and topological constraints, until the whole system is 
in the u = — 1 state. The energy of the initial state is 
given by E = (1/2) J dxp 2 , and the noise field thus has to 
be assigned initially in order to reach the switched state 
u = — 1 in a prescribed time T. For topological reasons 
the domain walls must nucleate and annihilate in pairs in 
general accompanied by absorption or emission of linear 
modes. Since the linear modes also carry positive ac- 
tion the dynamical modes with lowest action correspond 
to nucleation or annihilation of domain wall pairs with 
equal and opposite momenta, i.e., equal speeds. 

In the case of periodic boundary conditions the mo- 
mentum IT = J dxudp/dx of the initial and final states 
is zero. The system is translational invariant and the 
formation and annihilation of one or several domain wall 
pairs moving with the same speed take place at equidis- 
tant positions along the axis. For fixed boundary con- 
ditions the translational invariance is broken and the 
momentum IT is nonvanishing corresponding to nucle- 
ation and annihilation of domain walls at the boundaries. 
This general scenario of switching is completely consis- 
tent with the numerical analysis in Q. 

Switching a system of size L in time T by means of a 
single domain wall, corresponding to the pathway via the 
lowest local saddle point of the free energy at i^w = m, 
the propagation velocity v — po/m — L/T and we obtain 
the action S±(T) — mL 2 /2T and associated transition 
probability P oc exp(— mL 2 /2AT). In the thermody- 
namic limit L— >oo,P^0asa result of the broken 
symmetry in the double well potential. At long times 
the action falls off as 1/T. At intermediate times t and 
positions x we have P oc exp(— mx 2 /2At) and we in- 
fer that the domain walls in the stochastic interpretation 
perform a random walk with mean square displacement 
2Ai/m, corresponding to diffusive behavior. 

In Fig. |l|a we show a domain wall nucleating at the 
left boundary and propagating with constant velocity v = 
1/T to the right boundary, where it annihilates. We have 
used the same parameter values as in Q , i.e., 5 = T = .03, 
2Tk 2 = S" 1 , T = 7, and a system size L = 1. In Fig. [l]b 
we have plotted the trajectory of the domain wall in space 
and time. The switching can also take place by nucleating 
two domain walls at the boundaries. These then move at 
half the velocity v/2 and subsequently annihilate at the 
center. This process corresponds to the pathway via the 
local saddle point of the free energy at Fd w = 2m, and the 
action is given by S%(T) = 2Si(4T). Snapshots of this 
process are shown in Fig. ^|a and the corresponding space- 
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FIG. 1: In a) the switching from u = +ltou = — lin time 
T is effectuated by means of a domain wall propagating with 
velocity v — 1/T. The domain wall is nucleated at x = and 
annihilated at x = 1. In b) the process is depicted in an (x, t) 
plot. 



time plot in Fig. Combining the contributions from 
nucleation and the subsequent domain wall propagation, 
we can write a heuristic expression for the total action: 



S n {T) = nS n 



mL 2 
2nT' 



(12) 



where S nuc is t ne action for nucleating a single domain 
wall and n is the number of walls. The action of nucle- 
ation is easily estimated from the Arrhenius factor asso- 
ciated with the Kramers escape from the ground state 
u = +1 to the saddle point in the free energy (|J), i.e., 
P oc exp(— rFdw/A). We thus obtain a nucleation ac- 
tion 5 nuc of order Tko. A more detailed argument, to be 
presented elsewhere, based on estimating S nuc from the 
equations of motion @ and (0) and the action (0) , yields 
for domain wall nucleation and annihilation 



Sn 



6.5rfc . 



(13) 



In Fig. |^ we have plotted S versus T for n = 1 — 6 domain 
walls using the parameter values in [Q. Choosing S^uc 
according to (|l3|) we find excellent agreement with the 
numerical results. As also discussed in j| we note that 
the switching scenario depends on T. At shorter switch- 
ing times it becomes more favorable to nucleate more 
domain walls. In the present formulation this feature is 
associated with the finite nucleation or annihilation ac- 
tion iSnuc. This is evidently a finite size effect in the 
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FIG. 2: In a) the switching from u = +1 to u — — lin time 
T takes place by means of two domain walls propagating in 
opposite directions with velocity v = 1/2T. The domain walls 
nucleate at the boundaries and annihilate at the center. In 
b) the switching process is depicted in an (x,t) plot. 




FIG. 3: The action S(T) given by (|12|) is plotted as a function 
of T for transition pathways involving up to n = 6 domain 
walls. The lowest action and thus the most probable transi- 
tion is associated with an increasing number of domain walls 
at shorter times, indicated by the heavy limiting curve. The 
curves correspond to choosing S n uc = 5rfco. 



sense that the action at a fixed T diverges in the ther- 
modynamic limit L — > oo, corresponding to the broken 
symmetry. 



In this letter we have presented a dynamical descrip- 
tion and analysis of a specific non equilibrium transition 
in the noisy Ginzburg-Landau equation based on a canon- 
ical phase space formulation. We find good agreement 
both qualitatively and quantitatively with the numerical 
finding of E. et al. Q] based on an optimization of the 
Frcidlin-Wentzel action. The dynamical approach offers 
in the nonperturbative weak noise or low temperature 
limit an alternative way of determining dynamical path- 
ways and the Arrhenius part of the associated transition 
rates. 
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